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(N| ■ ABSTRACT 



Context. The Fundamental Plane relates the structural properties of early-type galaxies such as its surface brightness and effective 
radius with its dynamics. The study of its evolution has therefore important implications for models of galaxy formation and evolution. 
Aims. This work aims to identify signs of evolution of early-type galaxies through the study of parameter correlations such as the 
Fundamental Plane, the Kormendy, and the Faber-Jackson relations, using a sample of 135 field galaxies extracted from the Extended 
£N| 1 Groth Strip in the redshift range 0.2<z<1.2. 

Methods. Using DEEP2 data, we calculate the internal velocity dispersions by extracting the stellar kinematics from absorption line 
spectra, using a maximum penalized likelihood approach. Morphology was determined through visual classification using the V+I 
images of ACS. The structural parameters of these galaxies were obtained by fitting de Vaucouleurs stellar profiles to the ACS I-band 
images, using the GALFIT code. To check the effect on the Fundamental Plane of the structural parameters, Sersic and bulge-to-disc 
decomposition models were fitted to our sample of galaxies. A good agreement was found in the Fundamental Plane derived from the 
. three models. 

P> l Results. Assuming that effective radii and velocity dispersions do not evolve with redshift, we have found a brightening of 0.68 mag 

in the B-band and 0.52 mag in the g-band at <z>=0.7. However, the scatter in the FP for our high-redshift sample is reduced by half 
when we allow the FP slope to evolve, suggesting a different evolution of early-type galaxies according to their intrinsic properties, 
such as total mass, size or luminosity. The study of the Kormendy relation shows the existence of a population of very compact (R e <2 
Kpc) and bright galaxies (-21.5>M^>-22.5), of which there are only a small fraction (0.4%) at z = 0. Studying the luminosity-size 
and stellar mass-size relations, we show that the evolution of these compact objects is mainly caused by an increase in size that 
could be explained by the action of dry minor mergers. We detect also an evolution in the Fundamental Plane caused mainly by this 
population of very compact and bright galaxies. Unfortunately, we cannot distinguish a change in the slope from an increase in the 
scatter of the Fundamental Plane since our high-redshift sample is biased to the brightest objects. 

Key words. Galaxies: evolution - Galaxies: fundamental parameters - Galaxies: elliptical and lenticular, cD - Galaxies: kinematics 
and dynamics 
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\ 1. Introduction ness and the effective radius, as defined by de Vaucouleurs, are 

closely correlated (the Kor men dy relation, KR). Su bsequently, 

Elliptical galaxies, mainly composed by old star populations and iDjorgovski & Davisl (Il987l) and lDressler et al.l (Il987l) found that 

a small amount of gas and dust, represent an advanced stage of these two parameters plus the central velocity dispersion are the 

the galaxy formation as well as merging proces ses, making th e axes of a coordinated system where the galaxies define a plane 

k^H study of their properties particularly interesting (lGerhardll2006l) . known as 'Fundamental Plane' (FP). The smaller scatter of the 

^ ■ These galaxies are characterized by a de Vaucouleurs stellar pro- Fundamental Plane provides an improvement factor of 2 over 

C3 file, and by an elliptical distribution of their components sup- the FJR as a distance estimator. Since the velocity dispersion is 

ported by random motions rather than the ordered rotation of preventing the gravitational collapse, it can be used for estimat- 

spiral galaxies. ing the galaxy mass by applying the virial theorem. Thus, the 

w i i -a . . . j , 1 study of fundamental relations and their change with redshift 

The correlation between the luminosity and the central line- J ., , . £ . * • i i 

£ • 1 . j • t i _ , t -i can provide valuable information on the dynamical masses and 

or-signt velocity dispersion known as the Faber-Jackson re- w r i at ^ rn i ± . 

i rr~r — ^7 — i inr^Tk i . ^ . n i . mass evolution of galaxies. Moreover, the FP relates structural 

lation (FJR; Faber & Jackson 1976) relates the stellar and to- c . . u . , , ... ^ u 

A , ' 7rr~. — : i • — . ^ ^ 1t properties of galaxies, such as sizes and luminosities with the 

tal mass of an elliptical galaxy m a similar way to the Tully- j u v 1 j vi j i i - 

, w . r • i j i u i "! p i i dynamics, hence its study with redshift has important lmplica- 

Fisher relation lor spirals, and has been used to rind relative / ^ ~ . \ w . r ,,. , , 

j . . 11- i i • / r t IT^ , Jhr>o^ k tions on the formation and evolution of elliptical galaxies. 

distances to elliptical galaxies (see lor example, Dressier 1984j). r ° 

However, while the properties of spirals are rather well described Local studies show that the FP is followed by both elliptical 

by the two quantities (luminosity and rotation velocity) involved (E) an d lenticular (SO) galaxies dJorgensen. Franx & Kjaergaardl 

in the Tully-Fisher relation, the large scatter in the FJR forced Il996l) . which are jointly referred to as early-type galaxies. At 

the authors to look for o t her pa rameters to describe the ellip- higher redshift, the FP of E and SO galaxies could be differ- 

tical galaxies. Kormend^ (1 19771) found that the surface bright- ent if the star formation activity of SO galaxies were more in- 
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tense in the past. However, some studies have compared the re- 
sults obtained for each type of galaxies and no differences be- 
tween elliptical and SO w ere found (e. g. Ivan Dokkum jfe Franxl 
U996l:lKelson et al. II 1997b . Nevertheless. lFritz et al.l (l2005l) found 
a mild and on average larger evolution of lenticular galaxies over 
ellipticals, investigating 96 galaxies in clusters at z~0.2. 

Pioneering evolutionary studies of the FP were made at in- 
termediate redshift (z~0.5) fo r ga laxies in clusters. Fo r instance, 
Ivan Dokkum & Franxl d!996b and lKelson etaD (1 1997b found that 
the FP at this redshift was similar to the FP for nearby galaxies, 
suggesting little changes on the structure of the older galaxies 
and a formation epoch at significantly higher redshift (zf orm >2). 
These authors also found a lower mass-to-light (M/L) ratio of 
the galaxies at z~0.5, consistent with passively evolving stel- 
lar populations, but the low number of galaxies in their sam- 
ples made difficult to quantify the results. In a subsequent work, 
Kel son etaD (l2000bl) increased the sample of galaxies in a clus- 
ter at z=0.33 to 53 galaxies, and derived mass-to-light ratios of 
the high redshift early-type galaxies lower than those in Coma 
cluster by A log(M/L v )=-0. 13+0.03. This change in the M/L v 
ratio would imply an increase of the V luminosity by -0.3 mag, 
if such evolution were due to th e passive evolution of st ellar 
populations alone. In this sense, van Dokkum et al.l (2001) and 
Zie gler et~al1 (120051) found a similar result from a larger sample 
of field galaxies, deriving an average increase of the B luminosi- 
ties by — 0.4 mag at <z>=0.4. In addition, [van Dokkum et al.l 
(2001) found that the M/L ratios of field galaxies evolves as 
Aln(M/L B )=(-1.35±0.35)z, similar to c luster galaxies. In con- 
trast, the studi e s at h igher redshift of iTreu et al.l ([2002) and 
iGebhardt et al.l ([2003) found a greater brightening for field 
galaxies than the one found for cluster galaxies. Restricting the 
IGebhardt et~ai1 (120031) sample to galaxies at z=0.4, the B lu- 
minosities of early-type galaxies would have been -0.7 mag 
brighter than nowadays, which was attributed to a more recent 
formation epoch, around z=1.5. 

While all previous au t hors found an evolution of the FP only 
in the offset, iTreu et all (120051) find a change in the slope of 
the FP which can be interpreted as a mass-dependent evolution. 
Then, massive galaxies evolve passively on longer time-scales, 
whereas less massive systems have more extended star formation 
histories and continue to form stars at much l ater epochs. The 
recent studies of iFritz. Bohm & Ziegled (120091) . for a sample of 
24 fi eld galaxies covering redshifts 0.20<z<0.75, and lFritz et al.l 
(2009), for a sample of 50 cluster galaxies at z~l, also support 
the evolution in the slope of the FP for both field and cluster 
galaxies, suggesting that the internal properties of a galaxy are 
more important to its evolutionary history than its enviroment. 

In the present work, we study the evolution of the FP of 
a sample of 122 early-type field galaxies at <z>=0.7 in the 
Extended Groth Strip (EGS). This larger sample of galaxies can 
provide more statistically significant results than other works, as 
well as to accurately determine the FP followed by high redshift 
galaxies. In addition, we analize a sample of 13 local galaxies 
following the same procedure used for the high redshift sample, 
to calibrate the local FP used as comparison. 

This paper is organized as follows. In Sect. 2, the data and 
sample selection criteria are presented. Determination of the ef- 
fective radius, suface brightness, and velocity dispersion are de- 
scribed in Sect. 3. In Sect. 4, the Fundamental Plane is derived 
and its evolution is analysed using the FJR and the KR. The last 
two sections provide the discussion and the summary of the re- 
sults. Throughout this article, the concordance cosmology with 
Qao = 0.7, D m o = 0.3 and Ho = 70 km s" 1 Mpc" 1 is assumed. 



Unless otherwise specified, all magnitudes are given in the AB 
system. 



2. Data & Sample Selection 

The sample consists of galaxies in the EGS sky region. The base- 
line for spectroscopy target pre-selection were the galaxies for 
which DEEP2 spectra (Data Release 3, DR3) in this field were 
available. We selected the objects, from the redshift catalogue, 
with the tags ZQUALITY=3 or 4 and CLASS =GALAXY to re- 
move stars and AGNs. The DEEP2 project dDavis et al.ll2003L 
2007) is a survey using the DEIMOS multi-object spectrograph 
( Fab er et al.l 120031) in the Keck telescope, to study the distant 
Universe. The grating used was the 1200 1/mm one, covering a 
spectral range of 6500-9100 A with a dispersion of 0.33 A/px, 
equivalent to a resolution R=/l/A/l=4000. 

T he photometric d ata used here are part of AEGIS sur- 
vey dDavis et al.l 120071) . B, R, and I- band photometry wer e 
taken with the CFH12K mosaic camera (ICuillandre et al.ll200ll) . 
installed on the 3.6-meter Canada-France-Hawaii telescope 
(CFHT). These magnitudes are included i n the DEEP2 photo- 
metric catalogue (Data Release 1, DRh lCoil et al.1 120041) and 
the magnitude errors from sky noise and redshift are available 
too. The data in the V-band (F606W) were taken from HST 
catalogue, and obtained using the ACS camera. Although we 
adopted the I-band data from CFH12K, we used also the im- 
ages in the I-band (F814W) obtained with the ACS, for deriving 
the structural properties of our sample of galaxies. In addition, 
we used the z-band magnitudes of the Canada-France-Hawaii 
telescope legacy survey (CFHTLS), since z-band at z=l roughly 
matches the rest-frame B-band photometry. B, R, and I-band 
mag nitudes are al r eady corrected for Galactic reddening based 
on ISchlegel et al.l (1 19981) dust maps. The V and z-bands have 
been corrected following the prescriptions of the same work. 
At this point, our pre-selected sample has 3862 galaxies with 
DEEP2 spectrum and B, V, R, I, and z-band magnitudes avail- 
able. 

The second selection cr iteria was the morphology. In 
iFernandez Lorenzo et al.l ([2009) we compared the DEEP2 spec- 
tra with a synthetic, early-type galaxy spectrum, with the aim 
of segregating E/S0 from non-E/S0,through the rms of the dif- 
ference between spectrum and template. We checked the result 
through a visual classification, and we found that the mismatch 
for the non-E/SO group was -3%, whereas the one associated 
with the group of E/S0 galaxies was significantly worse, due to 
contamination by spirals. Then, although we cannot select our 
sample of E/S0 this way, a preliminary selection using this pro- 
cedure allowed us to eliminate a large fraction of the sample 
(2/3), losing only a 3% of the E/S0 galaxies. For the remaining 
1303 galaxies, which are characterized by having a spectra sim- 
ilar to that of the template, we made a visual classification using 
the summed V+I images of HST/ACS, and we found 400 E/S0 
candidates. 

Finally, we selected the E/S0 galaxies that had, at least, one 
absorption line in the spectrum for determining the velocity dis- 
persion, such as E-band region (Fe line, AA5210), G-band re- 
gion (Fe and Ca lines, AA4300), and the H + K region (double 
Call line, AA3934, 3969). The double absorption line of sodium 
or the triple absorption line of magnesium were not used since 
they could not be deblended, due to the lack of spectral resolu- 
tion. 
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Fig. 1. Examples of HST/ACS surface brightness modelling. For each galaxy, we have represented two model fits: a single de 
Vaucouleurs function (a) and a combination of r 1/4 + exponential disc component (b). For each galaxy and model we have shown 
the original data (left panel), the model fit (central panel) and the residual image (right panel). 



3. Data Analysis 

The final sample of galaxies with photometric information in the 
B, V, R, I and z-bands, in the redshift range 0.2<z<1.2 and clas- 
sified as E/SO, consists of 135 galaxies. The FP is constructed 
from three of the four fundamental parameters of the E/SO galax- 
ies: luminosity, velocity dispersion, effective radius and surface 
brightness at the effective radius. Physical parameters of the 
galaxies in our sample are shown in table 1 (in electronic for- 
mat). Here, we describe the way to determine these parameters. 



3.1. Structural parameters 

Historically, the surface brightness and the effective radius 
of the early-type galaxies have been calculated by fitting 
the mean radial profiles with a de Vaucoleu rs function r 1/4 
(iDjorgovski & Davislll987l:lDressler et al.| [T987). In more recent 
evolutionary studies of the FP, other surface brightness profiles 
such as Sersic function, expon ential law, or combinatio ns of two 
types of models (for example. Ivan Dokkum et al .2001, find that 
bulge+disc models provide better fits than a pure r 1/4 law) have 
been used. Although the surface brightness and the effective ra- 
dius are strongly dependent on the fitting profile, different fits 
have little effect on the FP, since the galaxies shift a long the 
plane instead of perpendicular to it dKelson et al.ll2000bh . 

In the present work, the structural modeling of the galaxies 
were made on the I-band HS T/ACS images, using the GALFIT 
package of iPeng et al.l (120021) . The effects of colour gradients 



are minimised by using the measurement in the I-band, since it 
corresponds to an optical rest-frame band in our whole redshift 
range. GALFIT minimises de^ 2 residual between the data image 
and the model, adjusting free parameters simultaneously, such as 
the total magnitude or the effective radius. Since different fitting 
functions are available in this code, and it is possible to deblend 
objects into one or more components, we explored the results 
of three profile models: a pure de Vaucouleurs function, a Sersic 
profile, and a combination of r 1/4 and exponential laws. All mod- 
els were convolv ed with a p oint-spread function (PSF) gener- 
ated by Tinytim (Krist 1995) , and we used th e parameters de- 
rived by SExtractor dBertin & Arnoutsll 19961) in the combined 
V+I HST images as inputs in the GALFIT code. Fig. 1 presents 
examples of the results of the surface brightness modelling using 
a single de Vaucouleurs function and a bulge-to-disc decompo- 
sition (B/D). The mean difference in surface brightness between 
both fits is -0.3 mag. However, the FPs obtained from the struc- 
tural parameters calculated by t hese two models are equivalent, 
according to the above result o f Kel son et al. (2000b). We have 
found also the same FP using the structural parameters obtained 
by fitting a Sersic profile. 

For a proper comparison of galaxies at different redshifts, 
measured surface brightnesses need to be converted to a com- 
mon rest-frame band. To calculate these k-corrections, we need 
to know the spectral energy dist ribution (SED) of the galaxy. In 
iFernandez Lorenzo et al.l (!20Tol) we found that the templa te ob- 
tained by the routine kcorrect dBlanton & R oweis 2007) from 
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Fig. 2. Example spectra (not flux-calibrated) of our sample of galaxies. The thin black line represent the rest-frame spectral data. 
The thick blue line is the galaxy model spectrum obtained from the convolution between the velocity profile and the stellar templates. 
The dashed vertical lines demarcate the regions that were excluded from the fitting because of the sky contamination by molecular 
bands. 
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a non-negative linear com bination of five templates based on the 
iBruzual & Charlotl (120031) stellar evolution synthesis codes, pro- 
vides the best result for retrieving the SED of the galaxy. In ad- 
dition, we found that the most reliable k-correction is obtained 
from information in an observed band that roughly matches the 
rest-frame band. Here, we have the magnitudes in the B, V, R, I, 
and z-bands. Since the z-band roughly matches the B-band at 
z= 1 , these five bands provide enough photometric information 
to calculate the rest-frame B-band magnitudes in our redshift 
range. Since different instruments have been used to obtain the 
magnitudes in different wavelengths, the same aperture might 
not ensure a consistent fraction of light in each band, due to the 
PSF, the seeing or the pixel scale, that act spreading more or less 
the object. Then, we compared the k-corrections derived from 
the B, V, R, I, and z magnitudes with those obtained using the 
aperture magnitudes (in a 1.5" radius) provided by CFHTLS in 
the g, r, i, and z-bands, and we found good agreement between 
both k-corrections. 

The GALFIT code provides the effective semimajor axis (aso) 
rather than the circular effective radius (r e ). Since local works 
use stru ctura l parameters defined in a circular annuli, we used 
r e =a5o y/b/a, where b/a is the axis ratio, for a proper comparison 
with other samples. The surface brig htness SBe=-2.51ogio£ e 
was determined using (Peng et al. 20Q2) 

F tot = 2nrl?: e e k nk- 2n r(2n)b/a , (1) 

where F tot is the total flux in the rest-frame band (i.e. the flux 
already k-corrected), H e is the surface brightness within the ef- 
fective radius r e , &= 1.9992^-0.3271, T is the gamma function, 
b/a is the axis ratio, and n is the Sersic index (n=4 for the clas- 
sical de Vaucouleurs profile). The surface brightness was also 
corrected for cosmological dimming multiplying F tot by (1+z) 4 . 

3.2. Velocity Dispersion 

The observable, luminous component of early-type galaxies is 
mostly made of stars. Therefore, their spectra can be described 
as a sum of individual stellar spectra redshifted according to 
their line-of-sight velocities. Assuming that the spectrum of 
all stars is given by a single template, then the galaxy spec- 
tra can be described as a convolution between the model tem- 
plate that represents the stellar spectrum and a broadening func- 
tion that describes the internal kinematics. Early methods for 
deriving the line-of-sight velocit y dispersions (LO SVD) used 
both the Fourier fitting technique (Fra nx et aI1 [l989), where the 
broadening function is recovere d from a deconvoluti on pro- 
cess, or the correlation method (iTonry & Davisl fl979l) . where 
the galaxy spectrum is cross-correlated against an essentially 
noiseless stellar template. These techniques are very sensitive to 
template mismatch and assume a Gaussian form for the broad- 
ening function, a simplifying assumption not adequate for the 
kinematics of elliptical galaxies that frecuently present multi- 
component stellar structures. These limitations force the re- 
searchers t o look for better ways to measure the stellar kine- 
matics (see Ide Bruyne et al.l|2003[ for an o verview of the vari- 
ous methods). For example JB ended (Il990h derived a new tech- 
nique (Fourier correlation quotient, FCQ), which is a 
hybrid method that provides a detailed analysis of the shape of 
the broadening func tion, and that has been used in recent stud - 
ies of the FP feiegler et alJl2005t iFritz. Bohm & Ziegledl2009h . 
Thanks to the increase in the computational speed, the LOSVD 
has begun to be ob tained in the pixel space (iRix & Wh ite 1992; 
Kels on et"aI1l2000al) . A treatment of the problem in pixel space 



has the advantage that gas emission lines or bad pixels can easily 
be eliminated, the effects of noise are more easily incorporated, 
and allows an easier error estimate. 

In this work, the velocity dispersion, <x, has been calculated 
using the Penalized Pix e l-Fitti ng method (pPXF) developed by 
ICappellari & Emsellem ( 2004), that works in the pixel space. 
This software extracts the stellar kinematics from absorption- 
line spectra of galaxies, using a maximum penalized likelihood 
approach. The code creates a model galaxy spectrum by con- 
volving a template spectrum by a parametrized LOSVD. The 
template spectrum is a combination of stellar templates with dif- 
ferent metallicities and spectral types, and the function describ- 
ing the LOSVD is based on the Gauss-Hermite series. The pa- 
rameters of the LOSVD are fitted simultaneously by minimizing 
the^ 2 , which measures the agreement between the model and 
the observed galaxy spectrum. The maximum penalized like- 
hood formalism, allows obtaining a solution that reproduces the 
details of the actual profile when the signal-to-noise (S/N) is 
high, while an adjustable penalty term is added to the;^ 2 to bias 
the solution towar ds a Gaussian shape when the signal-to-noise 
(S/N) is low (see ICappellari & Emsellem! 12004. for a complete 
explanation). The ability of the pPXF method to fit a large set of 
stellar templates together with the kinematics, allows eliminat- 
ing the template mismatch problem. 

The one-dimensional DEEP2 spectra used in this work, were 
extracted from the two-dimensional ones by using the routine 
extract Id. pro that is part of the DEIMOS s pec2d code . 
We use Home's optimal extraction algorithm (see iHornell 19861) 
along the locus of constant lambda. In the spatial direction, we 
used an aperture equal to the half-light radius determined in 
the previous section unless that value is lower than 3 (~PSF) 
or larger than 10 pixels (the pixel scale of spectra is 0.117 arc- 
sec/pixel). For these extreme cases we adopted 3 and 10 pixels 
respectively (6 galaxies have an effective radius slightly larger, 
but the difference in the velocity dispersion is within the errors). 
Since -43% of the sample have a minimum of 3 pixels, we made 
aperture corrections for these objects (see below). The resulting 
ID spectra were corrected to rest-frame wavelength before call- 
ing pPXF. S/N of spectra have values between 5 and 40, with 
a median value of 10. In addition, for minimizing the fitting 
time, we limited the wavelength region used in the analysis from 
lc-250A to lc+250A, where lc is the absortion line selected in 
each spectrum for determining the velocity dispersion. In some 
cases, when the absortion line was at the beginning or the end 
of the spectral range covered by DEEP2, a shorter wavelength 
region had to be used during the fit. 

To determine the velocity dispersion for our wide redshift 
range (0.2<z<1.2), we need to use different absorption lines, 
such as E-band region (Fe line, AA5210), G-band region (Fe 
and Ca lines, /L14300), and the H + K region (double Call line, 
/L13934, 3969). The H + K region is, a priori, not very suitable 
for measuring the velocity dispersion because the Call lines are 
intrinsically broad. Since the H + K regio n is the most accessible 
kinematic diagnostic in galaxies at z~ 1 , Kobu lnicky & Gebhardtl 
( 2000) investigated the velocity widths measured from the H + K 
region and from the HI line (21cm), for a sample of late-type 
galaxies. They found a good correlation between both measures, 
by using a maximum-penalized likelihood approach and a set 
of stellar templates of 1 2 Galactic A m ain-sequence through 
K giant stars. In addition, Gebhar dt et al.l (2003) obtained good 
agreement between the internal kinematic measured separately 
over the G and the H + K regions, using the same proce- 
dure for a sample of early-type galaxies. Nevertheless, choos- 
ing an appropriate set of stellar templates is of the utmost im- 
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Fig. 3. Edge-on projection of the Fundamental Plane in the B-band fitting the galaxies at z<0.35 (left) and fitting the data at z>0.35 
(right). The black points represent the galaxies with z<0.35, and the open red diamonds are the objects with z>0.35. The galaxies in 
each redshift range are fitted using the parameters a and b obtained for local galaxies from the literature. The zero-point (ZP) and 
the rms scatter (cr tot ) are shown in both cases. 
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Fig. 4. Edge-on projection of the Fundamental Plane in the B-band. The black points represent the galaxies with z<0.35, and the 
open red diamonds are the objects with z>0.35. The parameters a, b, and the zero-point were obtained by fitting the galaxies at 
z>0.35, through a non -linear weighted least squares fit. The green squares added in the right panel are the early-type galaxies of 
Gebhardt et al. ( 2003) with redshifts in the range 0.3<z<1.0. 



portance to obtain a reliable ki nematic from the H + K region 
( Ko bulnicky & Gebhar dt 2000). The stellar templates used in 
this work have been t aken from the Indo -U.S. Library of Coude 
Feed Stellar Spectra (iValdes et al.ll2004l) . which cover a spectral 
range from 3460 to 9464 A, at a resolution of ~1 A FWHM and 
at an original dispersion of 0.44 A pixel -1 . From this library, we 
have selected a subsample of 18 stars with different metallici- 
ties and spectral types ranging from A to K. For a proper com- 
parison with the DEEP2 spectra, each stellar template is con- 
volved with the quadratic difference between the DEEP2 and the 
Coude Feed Stellar Spectra instrumental resolution. Finally, we 
excluded some regions from the fit because of the contamina- 
tion by molecular bands of the sky. In Fig. 2 we have represented 
some examples of DEEP2 spectra, the best galaxy model fitted 
by pPXF, and the regions excluded from the fit in some cases. 



We use Monte Carlo simulations to estimate measurement 
errors for the LOSVD extraction method. With this aim, we 
added noise to the original data and repeated the full mea- 
surement process for 100 realizations. Finally, we calculated 
the standard deviation from the distribution of values given by 
the simulation and used them as velocity errors. We found that 
~90% of the objects in our sample had errors lower than 10%. 

Since we have used spectra extracted from physical sizes that 
not correspond to the effective radius in some cases, the velocity 
dispersion of these objec ts must be aperture corrected. We have 
followed the method of iMehlert et al.l (120031) . who obtained the 
aperture corrections by computing velocity dispersion gradients. 
The median value of the <x-gradient mea sured by these authors i s 
slightly steeper than the one proposed by lJorgensen et al.l (Il995l) . 
but consistent within the error. Moreover, the difference in the 
aperture correction obtained from both methods is lower than 1 
km/s for all objects in our sample. 
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4. Results 

4.1. The Fundamental Plane 

The FP relates the effective radius r e , the central velocity disper- 
sion <x, and the surface brightness within the effective radius SB e 
in the following form: 

log r e = a log cr + b S B e + c (2) 

To derive the evolution of the FP we have divided our sam- 
ple in two redshift ranges: galaxies at z<0.35, with a mean value 
of 0.27, which can be considered as local galaxies since they 
follow a relation that is similar to that found in nearby clus- 
ters (Kel son et al.l l2QQ0bb : and a high-redshift sample, which 
is made of galaxies at z>0.35, with a mean value of 0.68. 
The local sample is composed only of 13 galaxies. This num- 
ber of galaxies is insufficient for determining the local FP pa- 
rameters but it can be used to calibrate the local FP used as 
comparison to the high-redshift galaxies, since both samples 
were analysed following the same procedure. We a dopt then, for 
the loc al c oefficients, the same values u sed by iGebhardt et al.l 
(120031) and lFritz. Bohm & Zieglerl (|2009 ) in the rest-frame B- 
band: a =1.25 and b =0.32, and we calculated the zeropoint 
c for the local and high-redshift samples, by using a non-linear 
weighted least squares fit. These a and b values were obtaine d by 
these authors from the loc a l samp les of iFaber et al.l {1989) and 
Sagl ia. Bender & Dressier! (1 1993b respectively, using the same 
cosmology that we adopted in the present work. In Fig. 3, we 
present the edge-on projection of the FP in the local (left) and 
high-redshift samples (right), in addition to the zeropoint and 
total rms scatter obtained in each case. For our local sample, we 
obtained a zeropoin t of c = -9.09 3 , in good agreement with the 
previous result of lGebhardt et al.l (120031) (-9.062); and the re- 
sulting rms scatter cr tot = 0.113 is not sign ificantly higher than 
the local value (-0.1 dexi lFritz et al.ll2005l) . However, we found 
a different FP intercept for the high-redshift galaxi es, according 
to previous studies of F P evolution (for example Gebh ardt et al.l 
2003; Fritz et al. 2009). This difference in the zero-point is usu- 
ally interpreted as a difference in surface brightness caused by 
luminosity evolution, under the assumption that all early-type 
galaxies evolve in the same way, i.e. the coefficients a and b 
are independent of the redshift. From the change in the zero- 
point of the FP we found a brightening of 0.68 mag in the B- 
ban d for early-type galaxies at <z>=0.7. Using the expression 
of Ivan Dokkum & Franxl (1200 ll) for the evolution of a single- 
age stellar population (Locl/(t-t/ orm )^) and the same value of 
&s=0.91 used by these authors, this evolution of AM#=-0.68 
mag implies a formation redshift higher than z=5. This evolu- 
tion could be underestimated in the frame mark of progenitor 
bias effect. At z=0, we are observing early-type galaxies that 
could come from merged spirals at high-z, unobserved in high- 
redshift samples. Therefore, the high redshift sample is a bi- 
ased subset of the low redshift one. Ivan Dokkum & Franx (120011) 
studied the observed and true evolution of log(M/L B ) and found 
a difference of A log(M/L B )~0. 1 at <z>=0.7. Assuming that this 
difference is only due to luminosity, the true absolute magnitude 
is ~0.25 mag brighter than the observed one. Then, the luminos- 
ity evolution found in this work from the FP is underestimated. 
Nevertheless, the rms scatter obtained in the high-redshift sam- 
ple under this assumption is double than the local value. Then, 
either there is an evolution in the FP rms scatter, or the high- 
redshift galaxies follow a FP that is tilted with respect to the 
local one. 

To investigate the last possibility, we recalculated the FP co- 
efficients without any previous assumption. The result is pre- 
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Fig. 5. Faber-Jackson relation for the early-type galaxies in our 
sample. The black points represent the galaxies with z<0.35, and 
the open red diamonds are the objects wi th z>0.35. The dashed 
black line is the local relation used by [Fritz, Bohm & Zieg ler 
( 2009), and the dot-dashed red line is the fit to their distant sam- 
ple (0.2<z<0.75). The solid red line is the weighted fit to our 
high-redshift sample. 

sented in Fig. 4. The resulting plane is tilted with respect to 
the local one, and it has a rms scatter of cr tot =0.112, in good 
agreement with the local value. In addition, this high-redshift 
FP presents a rms scatter which is half of that obtained using 
the local a and b coefficients. Since our high-redshift sample 
consists of a large number of objects, this result is statistically 
significant. To check these hig h-redshift coefficients of the FP, 
we have represented the data of Gebh ardt et al.l (120031) (the red- 
shift range is similar to that used in the present work), together 
with ours, in Fig. 4 (right). We have recalibrated the data of 
IGebhardt et al.l d2003) to the AB system, for comparison. The 
objects with log(Re)<0.5 are consistent with the high-redshift 
FP, but the larger galaxies (log(R e )>0.5) have a lower zero-point 
than the si milar objects in our s ample. However, the rms scat- 
ter of the Geb hardt et al] (120031) data for the whole sample is 
also substantially reduced when they are fitted using a=0.53 and 
6=0.17. 

The change of the FP with redshift seems to indicate a dif- 
ferent evolution of early-type galaxies according to their intrisic 
properties, such as total mass, size or luminosity. In the next sec- 
tions we analise the FJR and the KR to clarify which of these 
properties is responsible for the change found in the present 
work. We also study the bias caused by the limited luminosity 
of the high-redshift sample, that can be affecting the results. 

4.2. The Faber-Jackson relation 

The FJR links the luminosity to the velocity dispersion of 
galaxies. Assuming that the velocity dispersion is related to 
the total mass of the galaxy, the FJR allows us to study the 
evolution in luminosity with the mass for a given galaxy. In 
Fig. 5, we represented the FJR in the B-band for all the galax- 
ies in our samp le. As comparison, we show the local rela- 
tion obtained by Fritz, Bohm & Ziegler (2009) from the Coma 
galaxies of ISaglia. Bender & Dressier] (1 19931) (determined us- 
ing the same cosmological assumptions that have been adopted 
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Table 1. Parameters of the FJ and Ko rmendy relations obtain ed by fitting DEEP2 data in the redshift ra nge z>0.35. In addition, 
we sho w the local relations derived by lFritz. Bohm & Ziegleri (l2009h from the Coma cluster galaxies of ISaglia. Bender & Dressier! 
(fl993h . 



Faber- Jackson relation Kormendy relation 

Redshift range intercept slope cr tota i intercept slope cr tota i 
Local relation of Fritz, Bohm & Ziegler (2009) T22 ^831 19^95 2l8 
z>0.35 DEEP2 galaxies -5.730 + 0.033 -6.897 + 0.014 0.819 18.553 + 0.033 4.142 + 0.073 0.586 



in the present work), as well a s the fit to the distant sample 
of iFritz. Bohm & Ziegleri (120091) (0.2<z<0.75). We have recal- 
ibrated our data to the Vega system for comparison with the 
work of IFritz. Bohm & Ziegleri (l2Q09h . The galaxies in our lo- 
cal sample are in good agreement with the local relation, al- 
though our objects have in average lower velocity dispersions. 
However, for the high-redshift sample, we found that the galax- 
ies are brighter than their local counterparts for a given velocity 
dispersion. In table 2, we present the parameters of the FJ rela- 
tion obtained by fitting DEEP2 data in the redshift range z>0.35 
and its total scatter, together with the local relation derived by 
IFritz. Bohm & Zie gler ( 2009). We found evidence of a different 
evolution of the galaxies as function of their masses. Then, al- 
though all the galaxies were brighter in the past, we obtained a 
stronger evolution for the lower mass galaxies. Nevertheless, the 
large dispersion of this relation makes difficult to obtain a reli- 
able conclusion on the s lope change of the FJR This result is 
similar to that found by IFritz. Bohm & Ziegleri (120091) . but the 
average evolution in the B-band magnitude is more noticeable 

in our high-redshift sample (AM B lmag), probably because 

our average redshift is larger than theirs. 



4.3. The Kormendy relation 

The Kormendy relation (KR) between the surface brightness and 
the effective radius and its change with the redshift, allows us to 
study the evolution in luminosity as function of the galaxy size. 

The KR in the rest-frame B-band for our local and 
high-redshift samples of galaxies is represented in Fig. 6. As 
in t he previous sections, we use d the local relation derived 
by Fritz, Bohm & Ziegler ( 2009) from the local sample of 
Sagl ia. Bender & Dresslerl(ll993l) . as comparison, showing good 
agreement with our local sample (as in the previous section, we 
have recalibrated our data to the Vega system). Although the 
weighte d fit to our local sample provi des a steeper relation than 
those of IFritz. Bohm & Ziegleri (|2009), this discrepancy may be 
attributed to a single galaxy with a very weak surface brightness 
(less than 22 mag). For the high-redshift sample, we found evo- 
lution in the surface brightness for a given effective radius, in 
the sense that galaxies were brighter in the past, with a change 
more noticeable for the smaller galaxies. In table 2, we also show 
the parameters of the KR obtained by fitting DEEP2 data in the 
re dshift range z>0.35, together with the local relation derived 
by lFritz. Bohm & Ziegleri (120091) . Similar resul ts have been re- 
cent ly found for higher red shift galaxies by iToft et al] (120091) 
and iDamjanov et al.l (120091) in the i and r-bands respectively. 
However, the variation in the slope of the KR is attributed by 
these authors to selection effects, because the sample is biased 
to the brightest galaxies. Since we have the same selection ef- 
fects, and our high-redshift galaxies have absolute magnitudes 
brigher than the local ones, we need a local sample covering a 
wide range of luminosity for comparison with our high-redshift 
galaxies. 
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Fig. 6. Kormendy relation for the early-type galaxies in our sam- 
ple. The black points represent the galaxies with z<0.35, and the 
open red diamonds are the object s with z>0.35. The dashed blu e 
line is the local relation used by Fritz, Bohm & Zieg EH ([2009), 
the thin black line is the weighted fit to our local sample, and the 
thick red line is the weighted fit to our high-redshift sample. 



The largest sample of local galaxi es is likely the Sloan 
Digital Sky Survey database (SDSS). iBernardi et all (l2003l) 
derived the FP relations in the g, r, i, and z-bands from a 
magnitude-limited sample of nearly 9000 early-type galaxies in 
the SDSS, covering a redshift range of 0.01<z<0.3. The struc- 
tural parameters of the SDSS sample were obtained by fitting 
a de Vaucouleurs model to the observed su rface brightness pro - 
file that accounts for the effects of seeing (IBernardi et al.ll2003l) . 
For comparison with our high-redshift sample, we used the sur- 
face brightnesses and the effective radii obtained by these au- 
thors only in the g-band. Then, we need to calculate the rest- 
frame g-band magnitudes for our sample of galaxies. To this 
purpose, we used the k-corrections derived from the aperture 
magnitudes (in a 1.5" radius) of CFHTLS in the g, r, i, and z- 
bands (see Sect. 3.1), Finally, we determined the surface bright- 
ness within the effective radius in the g-band using equation 1 . 
In Fig. 7 we present the comparison of the absolu te magnitudes 
distr ibution in the g-band for the SDSS sample of lBernardi et al.l 
( 2003) and for our sample of galaxies. The absolute magnitudes 
of the SDSS galaxies were calculated from the surface bright- 
ness and the effective radius, by using equation 1. Both dis- 
tributions are very similar, showing that the sample of galax- 
ies in EGS that have DEEP2 spectra, is appropriately selected 
as representative of the whole field, and then both samples are 
comparable. In the upper panel of Fig. 8, we represented the 
KR in the g-band for our sample of galaxies together with the 
Sloan data. Again, we found good agreement for the galaxies 
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Fig. 7. Absolute magnitudes distribution in the g-band of all 
DEEP2 galaxies in our sample (dashed red line) and of galaxies 
in the SDSS sample of iBernardi etafl (120031) (solid black line). 

in the local sample. Note that the IBernardi et al.l (120031) sample 
is magnitude limited and therefore the fainter objects in our lo- 
cal sample, that correspond to the galaxies with the highest sur- 
face brightness and the lowest effective radius, are rather scarce. 
On the other hand, the evolution found previously for the high- 
redshift sample is not so clear here. For a proper comparison, 
we limited our high-redshift sample and the SDSS one to the 
same luminosity range, -21.5>M^>-22.5, which roughly cor- 
responds to -21.0>Mg>-22.0. We represented the results in the 
bottom panel of Fig. 8. The largest objects in our high-redshift 
sample are in good agreement with the local ones, and there- 
fore no evolution in size or luminosity is found for these ob- 
jects. However, we have a population of objects with low ef- 
fective radius and high surface brightness that is not present in 
the local sample. We checked the profile fitting of these "small 
size" galaxies when we found this population, almost inexistent 
in the local universe. The fit seems good in all cases, with a 
very low^ 2 . The R e is between 4 and 10 pixels (semimajor axis 
between 5 and 15 pixels) and the values of b/a are distributed 
along the whole range (0.3-1). The objects are at least 2 times 
larger than the mean PSF of the ACS band used, and more or less 
half of them are ellipticals and the other half are lenticulars. We 
calculated the number of objects with R e <2 Kpc and luminosi- 
ties -21.5>Mg>-22.5 that exist in the comoving volu me corre- 
spon ding to each sample. In the SDSS sample of IBernardi et al.l 
(120031) (V p =0.26 Gpc 3 ), we found -96 Obj/Gpc 3 , while the re- 
sult for our high-redshift sample (V p =0.001 Gpc 3 ) is -24000 
Obj/Gpc 3 , i.e. only the 0.4% of these objects with R e <2 Kpc and 
total luminosity -21.5>M^>-22.5 exist in the local universe. 
Then, an evolution in luminosity or size since z= 1 is necessary 
to explain the decrease in number of these bright and compact 
objects. 

5. Discussion 

In order to stud y the processes drivi ng early-type galaxy evo- 
lution at z< 1 , iMcIntosh et al.l (|2005) analysed the luminosity- 
size and stellar mass-size relations in the redshift range 0<z< 1 . 
They found evolution in the luminosity-size relation that is con- 
sistent with passive aging of the galaxy stellar population, with 
the largest evolution appearing for the smallest galaxies. In ad- 
dition, they ruled out any substantial evolution in the stellar 
mass-size relation for the same early-type galaxies, consistently 
with the passive evolution scenario. However, the evolution in 
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Fig. 8. Upper panel: Kormendy relation in the g-band for the 
early-type galaxies in our sample. The black points represent the 
galaxies with z<0.35, and the open red diamonds are the objects 
with z>0.35. The small yellow points represent the early-type 
galaxies of SDSS (IBernardi et al.l20 03). and the dashed blue line 
is the lineal fit to this whole sample. The thin black line repre- 
sents the weighted fit to our local sample, and the thick red line 
is the weighted fit to our high-redshift sample. In the bottom 
panel we have represented the same as for the upper panel, but 
limiting our high-redshift sample and the SDSS one to the ob- 
jects with absolute magnitudes in the range -21.5>M^>-22.5. 
The dot-dashed line correspond to M g =-21.5. 



the stellar mass -size relation has recently been stablish ed at 
higher redshifts (iTruiillo et al.l l2007h iBuitrago etaTll2008l) and 
the combination of Palomar massive galaxies sample and lower 
stellar masse s galax ies from the GEMS survey conducted by 
Truiillo et al. ( 2007) has shown a perceptible evolution in the 
stellar mass-size relation of early-type galaxies since z~0.65. 

In Fig. 9 we present the luminosity-size relation in the g- 
band (upper panels) and the stellar mass-size relation (bottom 
panels) for our sample of galaxies. The current surviving stellar 
mass of our galaxies is calculated by kcorrect as the model 
mass derived from the coefficients fit to each template. The 
stellar masses used in the present work are then assuming the 
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Fig. 9. Luminosity-size relation in the g-band (up) and stellar mass-size relation (bottom) for our sample of early-type galaxies. 
The black points represent the galaxies with z<0.35, and the open red dia monds are the objec ts with z>0.35, divided in two redshift 
ranges. The small yellow points are the early-type galaxies from SDSS dBernardi et al.ll2003l) and the dash ed line is the linear fit to 
these data. The solid and dashed lines in the bottom panels represent the stellar masses of ShenetaL|i2003), as function of the local 
size distribution of early-type galaxies and its 1 <x. 



Chabrier stella r initia l mass function. We used the SDSS data of 
iBernardi et al] (120031) as local comparison of the luminosity-size 
relation. For the local stellar mass-size relation, we overplotted 
th e distribution in th e stellar mass of SDSS early-type galaxies 
of IShen eta l. (2003) as function of the Sersic half-light radius. 
We found that the largest objects (R e larger than ~3 Kpc) follow 
very similar relations than their local counterparts. However, the 
smaller objects are on average brighter and have larger stellar 
masses than the z=0 objects with similar sizes. The same re- 
sult is obtained using the Sersic half-light radius determined by 
GALFIT for our sample of galaxies instead of the de Vaucouleurs 
one. Although a passive aging of the galaxy population, com- 
bined with different formation epochs, can explain the evolution 
in the luminosity-size relation, the change in the stellar mass 
is not compatible with this explanation. The evolution of these 
compact objects is then mainly driven by an increase in size. 
Ivan Dokkum et al] (120101) studied the growth of massive galax- 
ies since z=2 and found that it is due to a gradual buildup of 
their outer envelopes. These authors found that the mass in the 
central regions is roughly constant with redshift, whereas the 
growth of the outerparts is due to a combination of star forma- 
tion and mergers. Also, that star formation is important at the 
highest redshifts, but the galaxy growth is domitated by mergers 
at z<l, which is expected in ACDM galaxy formation models. 



Gas-rich mergers increase the star formation activity and can- 
not explain the old stellar ages of early-type galaxies. However, 
dry (dissipationless) mergers are more efficient increasing the 
size than the stellar mass, and therefore would be the domi- 
nant mechanisms of size and stella r mass growing. In addition, 
iNaab. Johansson & Ostrikerl (l2009h performed a hydrodynamic 
simulation and showed that the mass assembly histories domi- 
nated by minor ("dry") mergers and accretion of stars can ex- 
plain the increase in size whereas major mergers cannot. 

On the other hand, a d ecrease in velocity dispe rsion 
with redshift is obtaine d by Ivan Dokkum et al] d2010h and 
ICenarro & Trujillo ( 2009), m contrast to our apparent increase 
(see Fig.5). However, no evolution is found in the F JR if we com- 
pare o ur sample of galaxies and the local data of IBernardi et al] 
( 2003|) due to the large scatter of the local relation. Since our 
velocity dispersion is measured only in the central part of the 
galaxy, and this region is roughly constant with redshift, our re- 
sult of no change in the velocity dispersion is consistent with 
this scenario of "dry" minor mergers that act increasing the size 
of the galaxies, due to a buildup of their outer envelopes. In ad- 
dition, for galaxies with similar masses and luminosities, this 
change in size would be more noticeable for the most compact 
and small objects, while the growth in size for the largest objects 
would remain undetected. 
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Fig. 10. Edge-on projection of the Fundamental Plane in the g-band for the early-type galaxies in our sample. The black points 
represent the galaxies with z<0.35, and the open red diamonds are the objects with z>0.35. The small yellow points are the SDSS 
early-type galaxies of lBernardi et al] (120031) . represented using the FP coefficients obtained from the maximum likelihood ortogonal 
fit by these authors. In the left pa nel, we show the zer o-point (ZP) and the rms scatter (cr tot ) for our local sample when the a and b 
coefficients are fixed to those of lBernardi et al] (120031) . In the right panel we present the same as for the left panel, but limiting our 
high-redshift sample and the SDSS one to the objects with absolute magnitudes in the range -21.5>M^> -22.5. 



5.1. Effects of compact objects on the Fundamental Plane 

To study the effects on the FP of these population o f compact and 
brigh t objects, we used again the SDSS sample of lBernardi et al.l 
(2003). In the left panel of Fig. 10, we present the FP in the g- 
band of our data compared to the SDSS maximum likelihood 
ortogonal fit. We found good agreement for our local sample 
of galaxies and for the high-redshift sample, the result is very 
similar to that was found for the B-band. In this case, from the 
change in the zeropoint of the FP, we found a brightening of 
0.52 mag in the g-band for early-type galaxies at <z>=0.7. In 
the right panel of Fig. 10 we represent the same FP but restricting 
our high-redshift sample and the SDSS one to the objects with 
absolute magnitudes in the range -21.5>M^>-22.5. As in the 
rest of relations studied in the present work, the largest galaxies 
in this sample are distributed following the same FP as their local 
counterparts. Nevertheless, the smallest objects occupy an area 
of the plane that is poorly populated by local objects with the 
same characteristics. Note that the few objects of SDSS that have 
the same properties as our high-redshift compact objects, are 
distributed in a similar way in the FP, pointing to a decrease in 
this type of objects rather than other possible changes, such as 
luminosity differences caused by passive evolution. 

Although these compact and bright objects appear to be tilt- 
ing the FP, our high-redshift sample is biased to the brightest ob- 
jects. Therefore we cannot distinguish a change in the tilt from 
an increase in the scatter of the FP caused by these population of 
objects almost nonexistent in the local FP. 

6. Summary and conclusions 

We have investigated the evolution of the FP in the B-band, 
using a large sample of 135 early-type galaxies selected from 
the EGS in the redshift range 0.2<z<1.2. Morphology was de- 
termined through visual classification using the V+I images of 
ACS. The structural parameters of these galaxies were obtained 
by fitting de Vaucouleurs stellar profiles to the ACS I-band im- 



ages, using the GALFIT code. To check the effect on the FP of 
the structural parameters obtained using different profiles, Sersic 
and bulge-to-disc decomposition models were also fitted to our 
sample of galaxies. We found good agreement in t he FP derived 
from th e three models, a result already found by iKelson etaD 
( 2000b). Finally, the velocity dispersion was calculated by ex- 
tracting the stellar kinematics from the galaxy spectrum, using a 
maximum penalized likelihood approach. 

Assuming that effective radii and velocity dispersions do not 
evolve with redshift, we found a different FP intercept for the 
high-redshift galaxies that could be interpreted as a brighten- 
ing of 0.68 mag in the B-band and 0.52 mag in the g-band at 
<z>=0.7. However, the scatter in the FP for our high-redshift 
sample is substantially reduced when we allow the evolution 
of the FP slope, suggesting a different evolution of early-type 
galaxies according to their intrinsic properties, such as total 
mass, size or luminosity. 

To investigate the galaxy properties responsible for the evo- 
lution of the FP, we derived the FJ and the Kormendy relations 
for our sample of galaxies. We obtained a preliminary result of 
galaxies being brighter in the past for a given velocity dispersion, 
and evidence of a stronger evolution for galaxies with low ve- 
locity dispersions. However, the great dispersion of this relation, 
does not allow deriving reliable conclusions. For the Kormendy 
relation, we found evolution in the surface brightness for a given 
effective radius, in the sense that galaxies were brighter in the 
past, with a change more noticeable for the smaller galaxies. 
Since this evolution could be caused by selection effects, we 
compared a subsample of our high- redshift objects with a sub- 
sample of the SDSS local galaxies of lBernardi et al] (120031) cov- 
ering the same luminosity range. The surface brightness distri- 
bution as function of the effective radius of these objects is in 
agreement with those expected given their luminosity. However, 
we found a population of very compact (R e <2 Kpc) and bright 
galaxies (-21.5>M^>-22.5), of which only a small fraction 
(0.4%) exist at z = 0, and that is responsible of the apparent 
evolution in the Kormendy relation. Then, a change in luminos- 
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ity or size since z= 1 is necessary to account for the decrease in 
number of this bright and compact objects. 

For studying the processes driving the evolution of these ob- 
jects, we analysed the luminosity-size and stellar mass-size re- 
lations. We found that the smallest objects present higher lumi- 
nosities and stellar masses than their local counterparts, while 
the largest objects are similar to those at z=0. We ruled out 
passive evolution of stellar populations as responsible for this 
evolution since it cannot explain the change in the stellar mass. 
We think that the evolution of these compact objects is mainly 
caused by an increase in size, which could be explained by the 
action of "dry" minor mergers that act increasing the size of the 
galaxies more than the mass. 

Finally, we studied the effect of these bright and compact ob- 
jects in the FP by comparing the galaxies with -21 .5>M^>-22.5 
of our high-redshift sample and the SDSS one. The previ- 
ous evolution found in the FP seems to be caused mainly 
by these galaxies, which have virtually disappeared at z=0. 
Unfortunately, we cannot distinguish a change in the tilt from 
an increase in the scatter of the FP caused by these population of 
objects, since our high-redshift sample is biased to the brightest 
galaxies. 

Acknowledgements. This work was supported by the Spanish Plan Nacional de 
Astronomia y Astrofisica under grant AYA2008-06311-C02-01. We thank the 
DEEP2 group for making their catalogues and data publicly available. Funding 
for the DEEP2 survey has been provided by NSF grants AST95-09298, AST- 
0071048, AST-0071198, AST-0507428, and AST-0507483, as well as NASA 
LTSA grant NNG04GC89G. The work is based on observations obtained at the 
Canada-France-Hawaii Telescope (CFHT) which is operated by the National 
Research Council of Canada, the Institut National des Sciences de l'Univers of 
the Centre National de la Recherche Scientifique of France, and the University 
of Hawaii. This work is based in part on data products produced at the Canadian 
Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy 
Survey, a collaborative project of NRC and CNRS, and on observations obtained 
with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA. 

This study makes use of data from AEGIS, a multiwavelength sky survey 
conducted with the Chandra, GALEX, Hubble, Keck, CFHT, MMT, Subaru, 
Palomar, Spitzer, VLA, and other telescopes and supported in part by the NSF, 
NASA, and the STFC. 

The Millenium Simulation database used in this paper and the web applica- 
tion providing online access to them were constructed as part of the activities of 
the German Astrophysical Virtual Observatory. 

We thank the SAO/NASA Astrophysics Data System (ADS) that is always 
so useful. 



Fernandez Lorenzo, M., Cepa, J., Bongiovanni, A., et al. 2010, A&A, 521, 27 
Franx, M., Illingworth, G, & Heckman, T. 1989, ApJ, 344, 613 
Fritz, A., Ziegler, B. L., Bower, R. G, Smail, I., & Davies, R. L. 2005, MNRAS, 
358, 233 

Fritz, A., Bohm, A., & Ziegler, B. L. 2009, MNRAS, 393, 1467 

Fritz, A., J0rgensen, I., Schiavon, R. P., & Chiboucas, K. 2009, Astronomische 

Nachrichten, 330, 931 
Gebhardt, K., Faber, S. M., Koo, D. C, et al. 2003, ApJ, 597, 239 
Gerhard, O. 2006, Planetary Nebulae Beyond the Milky Way, 299 
Home, K. 1986, PASP, 98, 609 

Jorgensen, I., Franx, M., & Kjaergaard, P. 1995, MNRAS, 276, 1341 

Jorgensen, I., Franx, M., & Kjaergaard, P. 1996, MNRAS, 280, 167 

Kelson, D. D., van Dokkum, P. G, Franx, M., Illingworth, G. D., & Fabricant, 

D. 1997, ApJ, 478, L13 
Kelson, D. D., Illingworth, G. D., van Dokkum, P. G, & Franx, M. 2000a, ApJ, 

531, 159 

Kelson, D. D., Illingworth, G. D., van Dokkum, P. G, & Franx, M. 2000b, ApJ, 
531, 184 

Kobulnicky, H. A., & Gebhardt, K. 2000, AJ, 119, 1608 
Kormendy, J. 1977, ApJ, 218, 333 

Krist, J. 1995, Astronomical Data Analysis Software and Systems IV, 77, 349 
Mcintosh, D. H., et al. 2005, ApJ, 632, 191 

Mehlert, D., Thomas, D., Saglia, R. P., Bender, R., & Wegner, G. 2003, A&A, 
407, 423 

Naab, T., Johansson, P. H., & Ostriker, J. P. 2009, ApJ, 699, L178 

Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 

Rix, H.-W., & White, S. D. M. 1992, MNRAS, 254, 389 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Saglia, R. P., Bender, R., & Dressier, A. 1993, A&A, 279, 75 

Shen, S., Mo, H. J., White, S. D. M., et al. 2003, MNRAS, 343, 978 

Toft, S., Franx, M., van Dokkum, P., Forster Schreiber, N. M., Labbe, I., Wuyts, 

S., & Marchesini, D. 2009, ApJ, 705, 255 
Tonry, J., & Davis, M. 1979, AJ, 84, 1511 

Treu, T., Stiavelli, M., Casertano, S., M0ller, P., & Bertin, G. 2002, ApJ, 564, 
L13 

Treu, T., Ellis, R. S., Liao, T. X., et al. 2005, ApJ, 633, 174 
Trujillo, I., Conselice, C. J., Bundy, K., et al. 2007, MNRAS, 382, 109 
Valdes, F, Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 
251 

van Dokkum, P. G, & Franx, M. 1996, MNRAS, 281, 985 
van Dokkum, P. G, Franx, M., Kelson, D. D., & Illingworth, G. D. 2001, ApJ, 
553, L39 

van Dokkum, P. G, & Franx, M. 2001, ApJ, 553, 90 

van Dokkum, P. G, Pieter, G, Whitaker, K. E., et al. 2010, ApJ, 709, 1018 

Ziegler, B. L., Thomas, D., Bohm, A., et al. 2005, A&A, 433, 519 



References 

Bender, R. 1990, A&A, 229, 441 
Bernardi, M., et al. 2003, AJ, 125, 1866 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 

de Bruyne, V, Vauterin, P., de Rijcke, S., & Dejonghe, H. 2003, MNRAS, 339, 
215 

Bruzual, G, & Chariot, S. 2003, MNRAS, 344, 1000 (BC03) 

Buitrago, F, Trujillo, I., Conselice, C. J., et al. 2008, ApJ, 687, L61 

Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 

Cenarro, A. J., & Trujillo, I. 2009, ApJ, 696, L43 

Coil, A. L., Newman, J. A., Kaiser, N., et al. 2004, ApJ, 617, 765 

Cuillandre, J.-C, Starr, B., Isani, S., McDonald, J. S., & Luppino, G. 2001, 

ASPC, 232, 398 
Damjanov, I., et al. 2009, ApJ, 695, 101 

Davis, M., Faber, S. M., Newman, J., et al. 2003, SPIE, 4834, 161 
Davis, M., Guhathakurta, P., Konidaris, N., et al. 2007, ApJ, 660, LI 
Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 
Dressier, A. 1984, ApJ, 281, 512 

Dressier, A., Lynden-Bell, D., Burstein, D., et al. 1987, ApJ, 313, 42 
Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 

Faber, S. M., Wegner, G, Burstein, D., Davies, R. L., Dressier, A., Lynden-Bell, 

D., & Terlevich, R. J. 1989, ApJS, 69, 763 
Faber, S. M., Phillips, A. C, Kibrick, R. I., et al. 2003, SPIE, 4841, 1657 
Fernandez Lorenzo, M., Cepa, J., Bongiovanni, A., et al. 2009, A&A, 496, 389 



